Quantum theory of heating of a single trapped ion 
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The heating of trapped ions due to the interaction with a quantized environment is studied without 
performing the Born-Markov approximation. A generalized master equation local in time is derived 
and a novel theoretical approach to solve it analytically is proposed. Our master equation is in the 
Lindblad form with time dependent coefficients, thus allowing the simulation of the dynamics by 
means of the Monte Carlo Wave Function (MCWF) method. 
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Single ions confined in miniaturized radio frequency 
traps can be cooled down by means of laser cooling tech- 
niques very efficiently Q]. Experiments have shown that 
the center of mass of a single cooled trapped ion under- 
goes a quantized harmonic oscillatory motion Q . During 
the last decade a great deal of attention has been devoted 
to such systems since they turn out to be very weakly 
coupled to external environment, thus allowing the ma- 
nipulation and observation of the coherent dynamics of 
a single quantum system j^j. Indeed, by irradiating the 
ion with properly configured laser beams, it is possible 
to manipulate the vibronic state of the ion. The easi- 
ness to engineer at will the interaction between the mo- 
tional and internal degrees of freedom of the ion makes 
it possible to perform experimental tests of fundamental 
features of quantum mechanics. For example, many non- 
classical states of the oscillatory motion of the ion have 
been prepared and measured 0. By using multiple ions 
in a linear Paul trap, experiments on quantum nonlocal- 
ity have been performed [p | and many-particle entangled 
states have been realized [|. Furthermore cold trapped 
ions have been recently proposed as a physical implemen- 
tation for quantum computation |?| . 

Both for fundamental studies and for technological ap- 
plications it is of great interest to understand and study 
those factors limiting the fidelity of the operations per- 
formed to manipulate coherently the quantum state of 
the ion. Heating of the center of mass oscillatory motion 
of one or more trapped ions seems to be one of the ma- 
jor practical sources of decoherence [||. The process of 
heating of a single trapped ion is due to the electromag- 
netic coupling between noisy electric fields and the ion. 
Such fields give rise to fluctuating forces acting on the ion 
which can cause an increase in its motional energy. The 
physical origin of the noisy electric fields in the center 
of the trap has not yet been unambiguously identified. 
In fact, measurements of the heating rate are very dif- 
ficult to perform since they require high sensitivity and 
they may depend on parameters related to the specific 
trap geometry used || . For this reason a general theory 
of heating would help in identifying and possibly reduce 



the up to now unknown sources of noise. 

Most of the theoretical studies on the heating of 
trapped ions deal with the interaction between a single 
harmonic oscillator and a reservoir which is either de- 
scribed as a classical stochastic field H [To], [n], |l2| or 
as an infinite chain of quantum harmonic oscillators. In 
the first case an analytic solution for the dynamics of 
the system can be found without performing the Born- 
Markov approximation jnj. In other words, it is pos- 
sible to study the short time evolution of the heating 
function when the field noise spectrum is not flat, as it 
is indeed in any real experimental situation. However 
the analytical treatments, in this case, cannot describe 
the thermalization process, i.e. the long time behavior, 
since the stochastic field continuously feeds energy in the 
system, leading to an infinite growth of the heating func- 
tion. On the other hand, modelling a thermal reservoir 
as an infinite chain of harmonic oscillators at tempera- 
ture T allows, under certain approximations, to describe 
the dynamics of the system in terms of a master equation 
of the form O 



dp 

dt 



T(N + 1) [a)ap - 2apa) + pa) a 
TN \aa) p — 2a) pa + pacr] , 



(1) 



where a) and a are the creation and annihilation opera- 
tors of vibrational quanta, T is the heating constant and 
the environment is assumed to be at thermal equilib- 
rium, N being the mean number of reservoir excitations 
at temperature T. In Eq. ([!]) the terms proportional to 
T(N +1) account for spontaneous and stimulated emis- 
sion, that is a transfer of energy from the system to the 
reservoir, while the terms proportional to TN account for 
absorption, that is a transfer of energy from the reservoir 
to the system. The presence of the spontaneous emission 
term, arising from a quantized description of the reser- 
voir, ensures the thermalization process. Indeed, from 
equation (|l|) , it is possible to derive the following expres- 
sion describing the time evolution of the heating function: 
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(n(t)) = n(wo) (1 - e" rt ) , 



(2) 



with n(uj Q ) = (e 13 ^ - 1) , /3 = (i-TT) -1 and w fre- 
quency of the trap. In order to describe the dynamics of 
the system by means of the master equation (|l|), however, 
it is necessary to assume the validity of the Born-Markov 
approximation. As a consequence such a treatment is 
valid only for times t much bigger than the correlation 
time of the reservoir and thus does not give the correct 
short time behavior. 

The aim of our paper is to study the heating of a single 
trapped ion interacting with a quantized reservoir with- 
out performing the Born-Markov approximation. One of 
the main results of the paper is the derivation of an ana- 
lytical expression for the heating function describing both 
the non-Markovian short time behavior and the asymp- 
totic thermalization process. 

The Hamiltonian describing the interaction between 
the ion motion and a quantized reservoir modeled as an 
infinite chain of harmonic oscillators, can be written as 
follows H = Hy/a J^JLo c j (bj +bj) , where bj is the 

annihilation operator of reservoir excitations and a is the 
coupling strength. By using the time-convolutionless pro- 
jection operator technique Jj^l and performing the Rotat- 
ing Wave Approximation, we have derived the following 
master equation describing the non-Markovian dynamics 
of the system in the interaction picture: 



dp A(t)+ 7 (t) r _ . 

— — — [a)ap - 2apa) + pa) a] 

\aa) p — 2a) pa + paa^] . (3) 



dt 2 

A(t) - 7 (i) 



The functions A(t) and j(t) are defined in terms of the 
two reservoir characteristic functions |l5[ : the correla- 
tion function k(t) = a / °° cu\g(uj)\ 2 (2n(u}) + l) cos(uiT)dcu 
and the susceptibility p(r) = a oj\g(ui)\ 2 sin(u>T)dtu, 



where n(u)) 



Moreover, assuming a 



Lorentzian spectral density |g(^)| 2 = uj 2 /[tt(uj 2 + oj 2 )} 
(Ohmic environment |l6[ ) with tu c reservoir frequency 
cut, in the high temperature regime defined by the con- 
dition (h/3)^ 1 3> uj c and in the weak coupling limit, the 
real functions A(t) and j(t) can be written in a simple 
analytic form: 



A(t) 



-{e- Wot sin(a; c f/r) 



h(3ujQ r 2 
+ r[l-e- w "*coB(w c t/r3}, 



7 (t) = 



2 r 2 



l-e~~ Uat cos{uj c t/r) 



re 



- Wct wa.{uj c t/r) 



(4) 



(5) 



where r = lo c /loq. Remembering that in the high tem- 
perature regime n(w ) — (hfiuo) -1 , it is easy to check 



that in the Markovian limit, that is when t 3> tr with 
tr = l/u>c correlation time of the reservoir, equation (||) 
reduces to equation ([!]). The analytic expression of A(t) 
and 7(4) for a generic temperature T is, in general, more 
complicated. Nonetheless the temperature dependence of 
these coefficients of the master equation brings to light 
interesting features in the dynamics of the system and, 
for this reason, it is carefully discussed in a follow up 
paper 0. 

Our non-Markovian master equation (|^) has the gen- 
eral form of the time-convolutionless master equations 
]l8| , |l9| ] . However, it is worth noting that it is not only lo- 
cal in time but also, as long as the coefficients A(t)+j(t) 
and A(t) — j(t) are positive, in the Lindblad form. We 
have verified that, for the value of parameters currently 
used in the experiments and for r > 1 the time depen- 
dent sum and difference coefficients are positive at all 
times t. Note that the condition r > 1 simply means 
that the spectral density of the reservoir overlaps with 
the frequency of oscillation of the ion, which is a reason- 
able and commonly done assumption. In other words, we 
demonstrate for the first time that the non-Markovian 
dynamics of the reduced density matrix of the system 
here considered, is described by a Lindblad-type master 
equation. Thus, it is possible to use the standard MCWF 
method to unravel Eq. (||) |2(}]. We note that up to 
now, only quantum state diffusion unravelings have been 
known and used for simulating the temporal behavior of 
a harmonic oscillator interacting with a non-Markovian 
quantized environment f2l|| . We would like to emphasize 
that the advantage of working with a master equation of 
Lindblad type whose analytical solution may be success- 
fully found, as in our case, is that the positivity of the 
density matrix is ensured. In the following we will use 
the quantum jump MCWF method to obtain the tempo- 
ral behavior of the heating function and compare it with 
the analytic solution. 

We now briefly sketch the novel analytical treatment 
we have developed to solve the master equation (^|) and 
derive a closed expression for the heating function de- 
scribing both its short and its long time behavior. A de- 
tailed description of the method will be given elsewhere 
]l7t . The first step of our approach consists in expanding 
the density matrix of the system as follows |22] : 



P(t) = - / Xt(A, A*) exp (A*a - Aa f ) d 2 X (6) 



In Eq. (§), A € C and Xt(A,A*) = 

Tr {exp (Act^ — A*a) /5(i)} is the quantum charac- 
teristic function. Inserting Eq. (Q) into Eq. (|^) and 
by using some useful algebric properties of the crono- 
logically ordered time evolution superoperators T c (i) 
of the system, defined by p(t) — T c (t)p(0), we obtain 
the following expression for the quantum characteristic 
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FIG. 1: Short time behavior of the heating function for cou- 
pling constant a = 0.1Hz, trap frequency u>o = 10 7 Hz, r—W 
and T = 300K. We compare the analytic solution and the 
MCWF simulation with 10 r histories. 



140 




t(<») 



x 10 



FIG. 2: Long time behavior of the heating function showing 
the thermalization process, for coupling constant a — 0.1Hz, 
trap frequency ujo = W 7 Hz, r=10 and T = 10 mK. We com- 
pare the analytic solution and the MCWF simulation with 
10 4 histories. We note that, even for such a small value of the 
reservoir temperature, the high temperature approximation 
holds since [hpy 1 ~ 10 9 #2 which is one order of magnitude 
bigger than uo c . 



function: 

Xt (A,A*)=cxp(-A(i)|A| 2 ) X o(T 7 1 / 2 (t)A;r 7 1 / 2 (i)Aj. (7) 

In this equation xa is the quantum characteristic func- 
tion at the initial time instant t = 0, 



and 



T 7 (i) = exp (- 2 J^ l{t')dt^j 



A(t) = / T,{t)T-\t')K{t'). 
Jq 



(8) 



(9) 



The functions A(t) and j(t) are defined by eqs. (|J) and 
(||). For example, assuming as initial state a vibrational 
Fock state |fc), Eq. (Q) becomes 



Xt{\ A*)=exp 



A(t)- 



T 7 (i) 



|A| 2 



(10) 

where Lk (y) is the Laguerre polynomial of order k in y. 
Eq. (0) together with Eq. (^|) gives the analytic expres- 
sion of the density matrix of the system in an operatorial 
form at any time instant t . From Eq. (]?]) it is easy to 
derive the following expression for the mean vibrational 
quantum number, i.e. for the heating function: 

(n(t)) - A(i) + \ (T 7 (t) - 1) + fcT 7 (t). (11) 

In Figs. |l| and || we show the short and long time be- 
havior of the heating function of a trapped ion initially 
prepared in its vibrational ground state. We compare the 
analytical solution and the MCWF simulation performed 
starting from the master equation (^) . The figures show a 
very good agreement between the analytical and the nu- 
merical approaches. We note that our quantum theory 
of heating predicts the initial quadratic behavior typical 
of non-Mar kovian dynamics. In Fig. ^ we compare the 
non-Markovian and Markovian time evolution for times 
t smaller than the reservoir correlation time 7~r = 1/lo c . 
A similar result was deduced by James and by Budini 
in the case of a classical environment described in terms 
of a stochastic noisy electric field || O. Our approach, 
however, allows also to describe the asymptotic thermal- 
ization process as shown in Fig. ||. Indeed, for times 
t ^> tr and for a flat reservoir spectrum r> 1, the heat- 
ing function given by equation (11) reduces to the well 
known Markovian equation given by Eq. (||), as one can 
easily derive substituting the asymptotic expressions for 
A(t — > oo) and j(t — > oo) into equations (||), (51) and 
([fl]), in the limit r > 1. 

In conclusion we have proposed a new approach for 
studying the heating of a single trapped ion due to the in- 
teraction with a quantized reservoir. Our approach does 
not rely on the Born-Markov approximation and thus al- 
lows to describe the initial quadratic behavior of the heat- 
ing function. Moreover, due to the quantum description 
of the environment, our solution describes correctly the 
thermalization process. Stated another way our quantum 
theory of heating bridges a gap existing in the literature. 
Until now, indeed, only non-Markovian analytical solu- 
tions, in the case of interaction with a classical environ- 
ment, or Markovian analytical solutions, for interactions 
with a quantum environment, were known. 

It is important to emphasize that measurements of the 
initial quadratic behavior of the heating function are dif- 
ficult to perform for two reasons. Firstly because this 
would mean to make many measurements in an interval 
of time of the order of tr ~ 10 -8 s, for the value of pa- 
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FIG. 3: Comparison between the short time Markovian and 
non-Markovian time evolutions of the heating function for 
coupling constant a — 0.1Hz, trap frequency ujq — 10 7 Hz, 
r=lQ and T = 300-fsT. 



rameters used in the paper. The second and more impor- 
tant reason is related to the difficulty in revealing so small 
variations in (n(t)} as the ones shown in Fig. |l|. Very re- 
cently, however, it has been experimentally demonstrated 
the possibility of engineering both the type of reservoir 
interacting with a single trapped ion and the coupling 
between the system and the environment ^3|. The an- 
alytical solution we have presented in this paper makes 
it possible to look for ranges of the relevant parameters 
of both the system and the reservoir in correspondence 
of which the non-Markovian quadratic behavior becomes 
experimentally observable. 

The experimental ability in engineering reservoir sug- 
gests also another application. As we have already men- 
tioned in the paper, we have analyzed the dependence of 
the coefficients of our master equation (|J) on reservoir 
parameters such as its temperature T and frequency cut 
lu c . We have found that by changing such parameters 
the master equation passes from Lindblad type to non- 
Lindblad type. Such a modification does reflect a deep 
change in the dynamics of the system. This result will 
be discussed in more detail in a follow up paper since 
we believe that the possibility of studying analytically 
such changes may give more insight in understanding the 
fundamental properties of the heating process of single 
trapped ions. 

We acknowledge K.-A. Suominen for helpful comments 
and the Finnish Center for Scientific Computing (CSC) 
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